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SUMMARY 


An  analysis  is  made  of  the  gust  response  (including  Lending  moment) 
of  an  airplane  having  the  degrees  of  freedom  of  vertical  motion  and  wing 
bending  flexibility  and  basic  parameters  are  established.  A  convenient^ 
but  accurate,  numerical  solution  of  the  response  equations  is  developed 
which  is  very  well  suited  for  making  trend  studies.  An  example  treated 
shows  results  which  are  in  very  good  agreement  with  the  results  obtained 
by  a  more  precise  but  more  lengthy  method. 

The  method  of  determining  a  gust  causing  a  known  response  is  indi¬ 
cated  and  a  procedure  is  given  for  determining  the  response  of  an  air¬ 
plane  directly  from  the  known  response  of  another  airplane  by  eliminating 
the  common  gust  condition. 


INTRODUCTION 


In  the  design  of  aircraft  the  condition  of  gust  encounter  has  become 
critical  in  more  and  more  instances,  mainly  because  of  the  ever-increasing 
fli^t  speeds.  Aircraft  designers  have  therefore  placed  greater  emphasis 
on  obtaining  rational  methods  for  accurately  predicting  the  stresses  that 
develop.  As  a  result,  the  number  of  papers  dealing  with  the  prediction 
of  stresses  in  an  aircraft  traversing  a  vertical  gust  has  significantly 
increased.  (See,  for  example,  refs.  1  to  9")  Many  of  the  papers  have 
treated  the  airplane  as  a  rigid  body  and  in,  so  doing  have  dealt  with 
either  the  degree  of  freedom  of  vertical  motion  alone  (refs.  5  to  8)  or 
with  the  degrees  of  freedom  of  vertical  motion  and  pitch  (refs.  7  and  9). 


^This  paper  is  a  revision  and  extension  of  a  paper  entitled  "The 
Determination  of  the  Response  Due  to  Gusts  of  One  Aiipilane  From  the  Known 
Response  of  Another  Airplane"  published  as  TN  No.  Structures  kO,  British 
R.A.E.,  June  19^9^  which  was  completed  by  Mr.  Houbolt  during  a  temporary 
tour  of  duty  with  the  Royal  Aircraft  Establishment.  Since  the  present 
paper  is  complete  in  itself,  no  further  reference  to  the  earlier  paper 
is  necessary. 
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Of  greater  concern  in  the  consideration  of  gust  penetration,  however,  is 

the  influence  that  wing  flexibility  has  on  structural  response.  This 

concern  has  two  main  aspects:  (l)  that  including  wing  flexibility  may 
lead  to  the  calculation  of  higher  stresses  than  would  be  obtained  by 
rigid-body  treatment  of  the  problem  and  (2)  that  wing  flexibility  may 
introduce  some  error  when  the  airplane  is  used  as  an  instrument  for 
measuring  gust  intensity.  Thus,  there  are  also  many  recent  papers  which 

treat  the  airplane  as  an  elastic  body.  In  most  of  these  papers  the 

approach  used  involves  the  development  of  the  structural  response  in 
terms  of  the  natural  modes  of  vibration  of  the  airplane  (refs.  2  to  4). 
Others  used  a  mohe  unusual  approach,  as  for  example,  reference  1,  which 
deals  with  the  simultaneous  treatment  of  the  conditions  of  equilibrium 
between  aerodynamic  forces  aind  structural  deformation  at  a  number  of 
points  along  the  wing  span.  Whatever  the  approach,  however,  the  main 
disadvantage  of  these  elastic-body  analyses  is  that  they  are  not  very 
well  suited  for  making  trend  studies  without  excessive  computation  time. 

In  the  present  paper,  the  case  of  the  gust  penetration  of  an  air¬ 
plane  having  the  degrees  of  freedom  of  vertical  motion  and  wing  bending 
is  considered.  Wing  bending  was  chosen  because  designers  have  expressed 
greater  concern  about  the  influence  of  this  flexibility  on  gust  response 
than  they  have  about  other  types  of  flexibility.  The  paper  has  the 
objective  of  trying  to  establish  some  of  the  basic  parameters  that  are 
involved  when  wing  bending  flexibility  is  Included  and  of  developing  a 
method  of  solution  which  is  fairly  well  suited  for  trend  studies  without 
excessive  computation  time.  Such  a  procedure  would  be  useful  in  evalu¬ 
ation  studies  which  are  intended,  for  example,  to  evaluate  the  effect  of 
such  factors  as  forward  speed,  spanwise  mass  distribution,  gust  length, 
and  gust  shape. 

The  equations  for  response  (including  accelerations,  displacements, 
and  bending  moments)  are  derived  and  the  basic  parameters  outlined.  An 
easy  numerical  solution  for  the  response  which  is  readily  handled  either 
by  manual  or  machine  methods  is  then  given.  The  inverse  of  the  response 
problem  is  considered  briefly;  that  is,  the  method  of  determining  the 
gust  causing  a  given  response  is  indicated.  Finally,  on  this  basis,  a 
procedure  is  outlined  whereby  the  response  of  one  airplane  may  be  found 
directly  from  the  known  response  of  another  airplane  without  the 
necessity  of  establishing  the  gust  causing  the  known  response. 


SYMBOLS 


a 


a 


n 


slope  of  lift  curve 

deflection  coefficient  for  nth  mode,  function  of  time  alone 
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A 

b 

c 

Co 

E 

f(s) 

F 

g 

H 

I 


Ly- 


m 


Me, 

Mn 

An 

P 


s,a 


S 


aspect  ratio  of  wing 
span  of  wing 
chord  of  wing 
chord  of  wing  midspan 
Young’s  modulus  of  elasticity 

nondimensional  gust  force,  /  ^1^(3  -  a)da 

J  0  u 

external  applied  load  per  unit  span 
acceleration  due  to  gravity 
distance  to  gust  peak,  chords 
bending  moment  of  inertia 

nondimensional  bending-moment  factor  =  Kj  ^  pUYM^^ 

aerodynamic  lift  per  unit  span  of  wing  due  to  vertical  motion 
of  the  airplane 

aerodynamic  lift  per  unit  span  of  wing  due  to  gust 
mass  per  unit  span  of  wing 

net  incremental  bending  moment  at  wing  station  j 

moment  of  wing  area  about  spanwise  station  under  consideration 

generalized  mass  of  nth  mode 

incremental  number  of  g  acceleration 

load  intensity  per  unit  spanwise  length 

distance  traveled,  t,  half-chords 

^o 

wing  area 
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t,T  time,’  zero  at  beginning  of  gust  penetration 

u  vertical  velocity  of  gust 

U  maximum  vertical  velocity  of  gust 

V  forward  velocity  of  flight 

W  total  weight  of  airplane 

y  distance  along  wing  measured  from  airplane  center  line 

w  deflection  of  elastic  axis  of  wing,  positive  upward 

%  deflection  of  elastic  axis  in  nth  mode,  given  in  terms  of  a 

unit  tip  deflection 

response  coefficient  based  on  a^, 

Ucq 

a  second  derivative  of  Zq  with  respect  to  s 

P  second  derivative  of  Z]_  with  respect  to  s 

7^  bending-moment  response  factor,  ratio  of  bending  moment 

obtained  for  airplane  considered  flexible  to  bending  moment 
obtained  for  airplane  considered  rigid 

e  distance  interval^  half-chords 


nn 


P 


reduced- frequency  parameter. 


2V 


nondimensional  relative-density  parameter. 


apc^S 


nondimensional  bending-moment  parameter. 


apCoMcQ 


mass  density  of  air 

function  which  denotes  growth  of  lift  on  rigid  wing  entering 
a  sharp-edge  gust  (Kussner  function)  ^ 
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0^  natural  circular  frequency  of  vibration  of  nth  mode 

1  -  C  function  which  denotes  growth  of  lift  on  an  airfoil  following 
a  sudden  change  in  angle  of  attack  (Wagner  function) 

e  =  1  -  <I> 

Subscripts: 

j  spanwise  station 

n  natural  modes  of  vibration 

m  number  of  distance  intervals  traveled 

Notation: 


column  matrix 


square  matrix 


Dots  are  used  to  denote  derivatives  with  respect  to  time;  primes 
denote  derivatives  with  respect  to  s  or  a. 


ANALYSIS 

Equations  for  Structural  Response 


Equations  of  motion.-  Consider  an  airplane  flying  horizontally  into 
vertical  gusts,  and  suppose  that  it  is  desired  to  include  wing  bending 
flexibility  in  determining  the  stresses  induced  by  these  gusts.  The 
problem  is  actually  one  of  determining  the  response  of  an  elastic  wing 
subject  to  dynamic  forces.  For  dynamic  forces  of  intensity  F  per  unit 
length,  the  differential  equation  for  wing  bending  is,  if  structural 
damping  is  neglected. 


— -  El  — -  =  -mw  +  F 
Sy2  ay2 


(1) 


where  w  is  the  deflection  of  the  elastic  axis  referred  to  a  fixed 
reference  plane.  The  task  of  determining  the  deflection  that  results 
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from  the  applied  forces  F  may  be  handled  conveniently  by  expressing 
the  deflection  in  terms  of  the  natural  free- free  vibration  modes  of  the 
wing.  With  regard  to  the  flight  of  an  airplane  through  gusts,  exami- 
mation  of  a  number  of  acceleration  and  strain  records  that  have  been 
taken  in  normal  flight  with  several  different  aircraft  shows  that  the 
response  to  gusts  is  composed  primarily  of  a  rigid-body  vertical  trans¬ 
lation  and  fundamental-bending-mode  excitation  of  the  wing.  Thus,  the 
assumption  is  made  in  the  present  analysis  that  the  response  may  be  given 
with  fair  accuracy  by  considering  only  these  two  degrees  of  freedom. 

This  assumption  is  probably  invalid  when  the  airplane  is  flying  near  the 
flutter  speed,  for  then  a  large  amount  of  coupled  bending-torsion  dis¬ 
placement  may  occur.  (See  ref.  3» ) 

The  wing  deflection  is  thus  assumed  to  be  given  by  the  equation 


w 


ao  + 


aiWi 


(2) 


where  w-^  is  the  deflection  given  in  terms  of  a  unit  tip  deflection 
along  the  elastic  axis  of  the  wing  for  the  fundamental  mode,  and  aQ 
and  a-j^  are  functions  of  time  alone.  In  this  form  aQ  denotes  the 

free-body  vertical  displacement  of  the  airplane  (in  this  case  the  dis¬ 
placement  of  the  nodal  points)  and  a-^  is  the  part  of  the  wing-tip 

deflection  which  is  associated  with  the  fundamental  mode,  as  illustrated 
in  the  following  sketch: 


The  use  of  symmetrical  modes  implies  that  only  the  symmetrical  gust  is 
to  be  considered  hereinafter. 

Substitution  of  equation  (2)  into  equation  (l)  yields 


=  -m 


(sq  +  .  F 


(3) 
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From  the  following  relation  which  expresses  the  condition  for  natural 
fundamental-mode  vibration 


equation  (3)  may  be  written 


a^a^^w-j_  =  -ni(aQ  +  a-j^W]^)  +  F 


where  is  a  natural  circular  frequency  of  vibration  of  the  funda¬ 

mental  mode.  If  this  equation  is  integrated  over  the  wing  span  and  use 
is  made  of  the  following  known  orthogonality  condition  of  the  free-body 
and  fundamental  modes: 


mw-^  dy  =  0 


the  following  equation  results: 


Mono 


where  Mq  is  the  aliqjlane  mass.  Now^  if  equation  (4)  is  first  multi¬ 
plied  through  by  w^^  and  then  integrated  over  the  wing  span  and  use  is 
made  of  equation  (5)^  the  following  equation  is  obtained: 


Mi'di  +  a^^M^a-|_  =  /  Fw-]_  dy  (7) 

^b/2 

where  M]_  is  the  generalized  mass  for  the  fundamental  mode^  that  is, 

rV2 

M]_  =  /  mwQ_‘^dy.  Equations  (6)  and  (7)  represent,  respectively, 

^-b/2 

the  equations  for  free-body  motion  and  fundamental  wing  bending  and  can 
be  solved  if  the  forces  F  are  known. 
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For  the  present  case  of  the  aii^lane  flying  through  a  gust,  the 
force  F  is  composed  of  two  parts:  a  part  designated  by  Ly  due  to 
the  vertical  motion  of  the  airplane  (including  both  rigid-body  and 
bending  displacements)  and  a  part  Lg  resulting  directly  from  the  gust 
(this  latter  part  is  the  gust  force  which  would  develop  on  the  wing 
considered  rigid  and  restrained  against  vertical  motion).  These  two 
parts  are  defined  (see  refs.  1  and  3)  in  the  equation  of  F  as  follows: 


F  -  Ly  +  Lg  =  - 1  pcV 


w|l  -  (ll(t  -  ■’’^dT  +  ^  pcV  U\lf(t  -  T)dT 


where  1  -  ^(t)  is  a  function  (commonly  referred  to  as  the  Wagner 
function)  which  denotes  the  growth  of  lift  on  a  wing  following  a  sudden 
change  in  angle  of  attack  and  for  two-dimensional  incompressible  flow 
is  given  by  the  approximation 


r  -1  -0.09^t 

|l  -  0(ty  =  1  -  0.l65e  ""  -  0.335e 


•O.^t 


and  ^(t)  is  a  function  (commonly  referred  to  as  the  Kussner  function) 
which  denotes  the  growth  of  lift  on  a  rigid  wing  penetrating  a  sharp- 
edge  gust  and  for  two-dimensional  Incompressible  flow  is  given  by  the 
approximation 


f—  -0. 26^t 

L^(t^  =l-0.5e  - 


An  additional  term  which  involves  the  apparent  air  mass  should  be  included 
in  equation  (8);  this  mass  term  is  inertial  in  character  and  may  be 
included  with  the  structural  mass  (see  ref.  1)  although  it  is  usually 
small  in  comparison.  The  lift-curve  slope  a  may  be  chosen  so  as  to 
include  approximate  over-all  corrections  for  aspect  ratio  and  compressi¬ 
bility  effects. 

If  w  as  given  by  equation  (2)  is  substituted  into  equation  (8) 
and  the  resulting  equation  for  F  is  substituted  into  equations  (6) 
and  (7),  the  following  two  equations  are  obtained  for  the  case  of  a 
uniform  spanwise  gust: 


T)dT 


(11) 


2E 
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2Mi  ..  a>,^2Mn 

— —  an  +  — = - =  an 

apVS  apVS  ^ 


a,)  [1  - 


n 


$(t  - 


dT  + 


T)dT 


(12) 


■  where  (because  of  mode  symmetry) 


nh/2 

8  =  2/  c  dy 
^0 


Pb/2 

Si  =  2  /  cw^  dy  )> 


pb/2 

S2  =  2  /  cw,  ^dy 
^0 


(13) 


Equations  (ll)  and  (12)  may  be  put  in  convenient  nondimen sional 
form  by  introducing  the  notation 


or 


(14) 


-n  =  0^  =n  (15) 

where  Cq  is  the  midspan  chord  of  the  wing  and  U  is  the  maximum 
vertical  velocity  of  the  gust.  With  this  notation,  equations  (11) 
and  (12)  may  be  written 


+  rn  z 


1^1 


i|r(s  -  a)da 


(16) 


4k 
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+ 


da  + 


\lf(s  -  CT)(ia 


(17) 


where 


8M0 

^0  spCqS 


8M-, 

n  =  I.  — 

1  apCoS 


(18) 


and  a  prime  denotes  a  derivative  with  respect  to  a.  Equations  (16) 
and  (17)  are  the  basic  response  equations  in  the  present  analysis.  The 
five  parameters  appearing  in  these  equations  and  given  by  equations  (18) 
depend  upon  the  forward  velocity^  air  density,  lift-curve  slope,  and  the 
aiiT»lane  physical  characteristics:  the  wing  plan  form,  wing  bending 
stiffness,  and  wing  mass  distribution.  Experience  has  shown  that  vari¬ 
ations  in  the  physical  characteristics  cause  significant  variations  in 
the  first  three  of  the  five  parameters,  while  the  last  two  vary  only  to 
a  minor  extent.  The  first  three  are  therefore  the  most  basic  param¬ 
eters;  M.Q  is  a  relative-density  factor,  frequently  referred  to  as  a 

mass  parameter,  and  is  associated  with  vertical  free-body  motion  of  the 
airplane;  li^,  similar  to  |j.q,  is  the  mass  parameter  associated  with  the 

fundamental  mode;  and  \  by  its  nature  may  be  interpreted  as  a  reduced- 
frequency  parameter  similar  to  that  used  in  flutter  analysis. 
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It  is  significant  to  note  that^  if  any  one  of  the  three  quanti- 
ities  Zq^  and  u  appearing  in  equations  (16)  and  (17)  is  specified 

or  known^  the  other  two  may  be  determined,  Thus^  if  the  gust  is  known, 
the  response  may  be  determined,  or  conversely,  if  either  Zq  or  is 

known,  the  gust  may  be  determined.  A  useful  equation  relating  Zq  and 
may  be  found  by  combining  equations  (I6)  and  (I7)  so  as  to  eliminate 
the  integral  dealing  with  the  gust.  The  result  is  the  equation 


which  is  used  subsequently. 

It  may  also  be  of  interest  to  note  that  in  effect  defines  a 

frequently  used  acceleration  ratio.  From  equations  (12)  and  (11),  the 
rigid-body  component  of  the  vertical  acceleration  may  be  written 


^0 


M 


or,  when  expressed  in  terms  of  the  incremental  number  of  g’s. 


z:\n 


_  ao  _  4VU  ^  M 

s  C^g  0 


An  acceleration  factor  AQq  based  on  quasi- steady  flow  and  peak  gust 
velocity  is  now  introduced  according  to  the  definition 

AngW  =  I  pSv2  ^ 


The  ratio 


An 

Arts 


is  thus  found  to  be 


An 
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Where  the  gust  shape  is  represented  analytically  and  the  unsteady- 
lift  functions  are  taken  in  the  form  given  by  equations  (9)  and  (lO), 
solution  of  the  response  equations  may  be  made  by  the  Laplace  transform 
method^  but  such  a  solution  is  more  laborious  than  desired.  Therefore, 
a  numerical  procedure  which  permits  a  rather  rapid  solution  of  the  equa¬ 
tions  has  been  devised  and  is  presented  in  a  subsequent  section.  It  may 
be  well  to  mention,  however,  that  the  response  equations  are  suitable 
for  solution  by  some  of  the  analog  computing  machines. 

Bending  stresses.-  The  bending  moment  and,  hence,  the  bending 
stresses  that  develop  in  the  wing  due  to  the  gust  may  be  found  as 
follows:  The  right-hand  side  of  equation ^(l)  defines  the  loading  on 
the  wing;  suppose  this  loading  is  noted  by  p,  then 


p  =  -mw  +  F 


By  use  of  equations  (2)  and  (8),  and  the  notation  of  equations  (l4) 
and  (15) ^  this  equation  becomes 

p  =  «ni  +  ^l"^l)  "  apcVuT  (zq"  +  ^l"^l)  -  ^)(s  -  a)1da  + 

Cq  J  q  L  J 


u’>l/(s  -  a)da 


If  the  moment  of  this  loading  is  taken  about  a  given  wing  station, 

say  y^,  the  following  equation  for  incremental  bending  moment  at  that 
J 

station  would  result: 

nb/2 

=  j  p(y-yj)<iy 

=  -^(Wmo^o"  -  apYuJ  (Mc^Zq"  +  + 


-  pVMc^  /  u'\l;(s  -  cT)da 


0 


(20) 
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where  the  M's  hearing  double  subscripts  are  first  moments  defined  as 
follows : 


and  is  the  station  being  considered.  Division  through  of  equa¬ 

tion  (20)  by  the  quantity  |  pYUMc^  gives  the  following  equation  which 
is  considered  to  define  a  bending-moment  factor  Kj  at  wing  station  y. 


The  factor  ^  pYUMc^  may  be  regarded  as  the  maximum  aerodynamic  bending 

moment  that  would  be  developed  by  the  gust  under  conditions  of  quasi¬ 
steady  flow  and  with  the  wing  considered  rigid  and  restrained  against 
vertical  motion  at  the  root.  The  bending-moment  factor  Kj  may  thus 
be  seen  to  be  the  ratio  of  the  actual  dynamic  bending  moment  that  occurs 
to  this  quasi-steady  bending  moment  and  therefore  may  be  regarded  as  a 
response  or  an  alleviation  factor. 
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A  more  convenient  form  for  the  bending-moment  factor  may  be 
obtained  by  solving  equations  (16)  and  (17)  simultaneously  for  the 

quantities 

substituting  these  values  into  equation  (22).  With  these  operations  the 
following  equation  results: 


r 

L  -  (I)(s  -  ay]  da  and  f  zq"  H 

L  -  ®(s  -  a^ 

Jq  ^ 

Jo 

K,  = 


M 


L 


^  I  pVUMcq 


^r2r]_  -  T2 


iri  -  rg 


^0  -  \]^o  — 


ri  -  r3 


- 


ri-  -  rg 


ri  -  r 


^1  ^1 


^1  -  ^2 


where 


(23) 


Mr 


■3  Mr 


3% 


^0  apc^Mj 


0 


8M, 


’ll  =  - 

1  apCoMc 


0 


> 


(24) 


It  is  seen  that^  when  bending  moments  are  being  determined^  three  addi¬ 
tional  basic  parameters  (eqs.  (24))  appear.  The  similarity  of  t)q 

and  to  and  is  to  be  noted;  first  moments  of  masses  and 

areas  are  involved  rather  than  masses  and  areas. 

Reduction  to  rigid  case.-  It  may  be  of  interest  to  show  the  reduc¬ 
tion  of  the  response  equation  to  the  case  of  the  airplane  considered  as 
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a  rigid  body.  Thus^  if  Z]_  is  equated  to  zero  in  equation  (l6)^  the 
following  equation  for  rigid-body  response  is  obtained: 

(25) 

If  z-^'  is  set  equal  to  zero  in  equation  (22)  and  use  is  made  of 

equation  {2^),  the  following  equation  for  the  bending-moment  parameter 
for  the  rigid-body  case  is  obtained 

^  (^0  “  ^0)^0"  (26) 

where  Zq"  is  the  acceleration  of  the  airplane  considered  as  a  rigid 
body. 

Matrix  Solution  of  Response  Equations 

In  this  section  a  rather  simple  numerical  solution  of  the  response 
equations  (16)  and  (17)  is  presented.  The  procedure  is  readily  adapted 
to  either  manual  or  punch-card  machine  calculations. 

The  derivation  proceeds  on  the  basis  that  the  response  due  to  a 
given  gust  is  to  be  determined.  The  airplane^  just  before  gust  pene¬ 
tration^  is  considered  to  be  in  level  flight  and  hence  has  the  initial 
conditions  that  the  vertical  displacement  and  vertical  velocity  are  both 
zero.  These  conditions  mean  that  Zq,  Zj_,  Zq' ,  and  z^*  are  all  zero 

at  s  =  0.  The  gust  force  can  be  shown  to  start  from  zero  and,  there¬ 
fore,  the  additional  Initial  conditions  can  be  established  that  Zq” 

and  z-j^  are  also  zero  at  s  =  0.  By  the  numerical  procedure,  solution 

for  the  response  at  successive  values  of  s  of  increment  e  will  be 
made  and,  for  the  case  being  considered,  it  is  found  advantageous  to 
solve  directly  for  the  accelerations  rather  than  the  displacements. 

In  order  to  make  the  presentation  more  compact,  the  following 
notation  is  introduced: 


3  =  Zl" 

0  =  1  -  f 
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and 


f(s) 


^  \lr(s  -  a)dcT 


(27t) 


With  this  notation^  equation  (I6)  would  appear  simply  as 


^0°  "  +  ^”(s)  (28) 


In  accordance  with  numerical-evaluation  procedures^  the  interval 
between  0  and  s  is  divided  into  a  number  of  equal  stations  of 
interval  e  so  that  s  =  me.  The  product  of  +  r3_p)  and  0(s  -  a) 

is  assumed  formed  at  each  station  and,  with  the  use  of  the  trapezoidal 
method  for  determining  areas,  the  unsteady-lift  integral  in  equation  (28) 
may  be  written  in  terms  of  values  of  a  and  p  at  successive  stations 
as  follows,  where  the  mth  station  corresponds  to  the  value  s: 


+  •  •  •  +  9iS,.i  +  1  9o%)  + 


^^l(®m-l^l  ®m-2p2  + 


©iP 


m-1 


h  Vm) 


(29) 


in  which  0q,  Q-^,  .  .  .  are,  respectively,  the  values  of  the 

1-0  function  at  s=0,  s=c,  ...  ^aQ  and  pQ  do  not  appear 

because  of  the  initial  conditions).  With  this  equation,  equation  (28) 

may  be  written  at  various  values  of  s  or  at  successive  values  of  mj 
the  result,  for  example,  for  m  =  1  is 

PoQi  =  -e0Qaj_  -  er2_0Q3p  +  f3_ 

and  for  m  =  2, 

'^o“2  ""  "^(^®1°1  ®0°2)  "  ^^l(^^l*^l  ^0^2)  ^2 


and  f2  are  the  values  of  the  gust-force  integral  at  s  =  e  and  s  =  2e.  The 
thus  formed  may  he  combined  in  the  following  matrix  equation: 
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According  to  the  derivation  presented  in  the  appendix,  the  value  of  at  s  -  me  may  be 

approximated  in  terms  of  the  past-history  value  of  z-^"  by  the  following  equation: 
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The  square  matrix  C  is  seen  to  be  similar  to  the  other  square 
matrices  in  that  it  is  triangular  with  all  the  elements  in  one 
column  made  up  of  the  elements  in  the  previous  column  moved  down  one 
row. 


An  equation  in  | 3 |  alone  is  obtained 
this  equation  into  equation  (30)  to  yield 


by  substituting  |a(  from 


*  [B]]|p|  =  WJM  =  |f|  (33) 


which  is  the  basic  response  equation  relating  ,p  (that  is  to  the 

gust  force.  This  equation  represents  a  system  of  linear  simultaneous 
equations  where  the  order  of  the  matrix  is  arbitrary;  that  is,  the  equa¬ 
tions  may  be  written  up  to  any  desired  value  of  s  =  me.  The  solution 
for  response  can  therefore  be  carried  on  as  far  as  desired.  Fortunately, 
the  equations  are  of  such  a  nature  that  simultaneous  solution  is  not 


required.  As  mentioned,  each  of  the  matrices  [a]  ,  [bJ  ,  and  [c]  is 

triangular  with  all  elements  0  above  the  main  diagonal  and  with  all 
elements  on  the  main  diagonal  of  each  matrix  equal;  therefore,  the  main 

diagonal  elements  of  [d]]  will  also  all  have  the  same  value  and  the 
elements  above  this  diagonal  will  be  0.  If  each  element  on  the  main 
diagonal  of  |d]  is  denoted  by  d]L  [p£|  is  "the  matrix  D  with 

the  main  diagonal  elements  replaced  by  O’s,  then 


[D]  =  di[l]  +  [Di 

With  this  equation,  equation  (33)  niay  be  written 


(34) 
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Expanded,  this  equation  has  the  form 


It  can  be  seen  that  a  step-by-step  solution  for  the  successive  values 

of  p  may  now  be  made;  that  is,  Pj_  is  solved  for  first  then,  with  P3_ 

established,  P2  is  solved  for,  and  so  on  as  far  as  is  desired.  With 

the  value  of  |p|  thus  established,  solution  for  |a|  may  now  be  made 
directly  from  equations  (32).  Values  of  the  displacements  Zq  and 
may  be  obtained  directly  from  a  and  P;  z-^  may  be  obtained  from  equa¬ 
tion  (31);  and  Zq  may  be  obtained  from  this  same  equation  with  p 
replaced  by  a. 

Some  mention  should  be  made  with  regard  to  the  selection  of  the  time 

interval  6.  A  rough  guide  to  use  in  selecting  g  can  be  obtained  by 

considering  X.,  which  appears  as  the  characteristic  frequency  in  most 
response  calculations.  The  period  based  on  this  frequency  would  be 

p-f 

Tg  =  Experience  has  shown  that  a  time  interval  in  the  neighborhood 

of  1/12  of  this  period  yields  very  good  results  (in  general  less  than 
1  percent  error) ;  accordingly,  a  reasonable  guide  in  choosing  e  would 

be  the  equation  e  Some  convenient  value  near  that  given  by  this 

equation  should  be  satisfactory;  in  general,  it  will  be  found  that  e 
may  be  1  or  greater. 

The  procedure  thus  outlined  provides  a  rather  rapid  evaluation  of 
the  response  due  to  a  prescribed  gust.  With  the  response  thus  evaluated 
the  bending  moment  at  any  value  of  s  or  the  complete  time  history  of 
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bending  moment  may  be  found  by  application  of  equation  (23).  It  should 
be  evident  that,  if  response  values  for  either  Zq"  or  z-^'  are  known, 

the  gust  causing  this  response  can  be  found  by  suitable  manipulation  of 
equations  (30)  and-  (32).  Thus,  if  Zq"  is  known,  P  in  equations  (3OI') 

and  (32b)  may  be  eliminated  to  give  the  equation 

[m  -  j  -  ifi 

Direct  substitution  of  Zq"  in  this  equation  allows  | f |  to  be  deter¬ 
mined.  In  most  practical  cases  the  second  term  in  equation  (30b)  con¬ 
tributes  only  a  small  amount  and  may  be  dropped  with  little  resulting 
error  in  the  gust  force.  The  equation  for  |f  |  is  then  simply 


[A]|a|  =  ffl 


Determination  of  Response  of  One  Aiip)lane  From 
Known  Response  of  Another  Airplane 

In  general,  a  given  gust  condition  produces  different  responses 
either  for  two  different  airplanes  or  for  the  same  airplane  with  dif¬ 
ferent  loading  conditions  or  forward  velocity.  It  would  be  expected, 
however,  that  the  response  of  the  two  aiip)lanes  could  be  correlated 
through  the  common  gust  condition.  This  correlation  may  be  demonstrated 
quite  easily  by  means  of  the  equations  given  in  the  preceding  section. 

The  case  to  be  treated  is  as  follows:  The  time  history  of  bending 
moment  due  to  a  gust  is  assumed  to  have  been  measured  in  one  airplane 
and  it  is  desired  to  calculate  directly  from  this  time  history  what  the 
bending  moment  due  to  the  same  gust  would  have  been  in  another  airplane. 
Although  the  derivation  is  presented  in  terms  of  bending  moment,  a 
similar  derivation  could  be  made  in  terms  of  either  accelerations  or 
displacements. 

If  use  is  made  of  equation  (3I)  to  write  the  successive  values  of 
the  displacement  z-]_  in  terms  of  the  accelerations,  the  bending-moment 

factor,  equation  (23),  may  be  written  in  terms  of  the  accelerations  alone 
and  the  following  matrix  equation  for  K  may  be  formulated: 
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With  the  use  of  equations  (32)  this  equation  may  be  written 


|k|  =  [eJIpI 

where 

[Ej  =  ^[C]  .  e[l]  .he2[G] 


(36b) 


(37) 
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in  which  is  the  identity  matrix  and 


[G] 


Substitution  of  the  value  of  3  as  obtained  from  equations  (36)  into 
equation  (33)  gives  the  following  relation  between  the  gust  forces  and 
the  bending-moment  parameter: 

|f|  =  [D][E]“^|Kjn|  (38) 

The  gust- force  matrix  |f|  (see  eq.  (27b))  may  be  expressed  in 
terms  of  the  gust  velocity  by  the  following  process:  It  is  assumed  that 
the  initial  vertical  velocity  of  the  gust  is  zero  and  that  successive 
values  of  gust  velocity  of  increment  €  are  designated  by  U]_,  U2^ 

u^  .  •  .  •  First-order  difference  equations  are  used  to  approximate  the 
slope  of  the  gust  velocity,  so  that,  in  general. 


t  _  ^+1  ~  ^-1 


If  this  equation  is  used  and  the  integral  equation  (27b)  is  handled  by 
the  trapezoidal  integration  method  similar  to  that  used  for  equation  (29), 
the  gust  force  may  be  written  in  terms  of  the  successive  values  of  gust 
gradient  so  as  to  form  the  following  matrix  equation: 
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fl 

^2 

_  1 
"  2 

"3 

^4 

1 

^2 


u 

U2 

U 

^3 

u 

u 

(39) 


^3  -  ^2 

^4  -  ^1^2  ^3  -  'I'l  "^2  'I'l 


where  ''1^2^  1^3  ..  .  are  successive  values  of  the  \tf  function. 

Substitution  of  this  equation  into  equation  (38)  allows  for  the  solution 
of  j^j  in  terms  of  the  parameter  K  as 


-  W^&>][E]""|K,n| 


-ll 


(40) 


where  [\tf]  is  the  square  matrix  in  equation  (39).  Different  airplanes 

flying  through  the  same  gust  will  experience  the  same  vertical  gust 
velocity  for  equal  absolute  distances  of  gust  penetration^  that  is^ 


(Vt) 


airpleine  1 


=  (Vt) 


airplane  2 


From  equation  (14),  then,  the  following  conditions  must  prevail: 


=  (="=0)2 

=  (^'=0)2 


(41) 


where  the  subscripts  1  and  2  denote  airplane  1  and  2,  respectively. 
Satisfaction  of  the  latter  condition  insures  that  the  gust  velocity  as 
given  by  equation  (40)  would  be  the  same  for  the  two  aiiplanes  being 
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compared.  This  common  gust  condition  may  therefore  he  eliminated  to 
yield  the  result 

If  it  is  assumed  that  K  for  airplane  1  is  known,  then  K  for  air¬ 
plane  2  may  he  written 

IkIs  =  _[E][D]-l[>|r]  2  W-^[D][E]-l|Ki]i  (k2) 

where  again  the  time  interval  chosen  for  the  two  airplanes  must  satisfy 
relation  (4l).  Thus,  if  the  bending  moment  due  to  a  given  gust  sequence 
is  known  for  one  airplane,  the  bending  moment  that  would  develop  in 
another  airplane  encountering  the  same  gust  sequence  can  he  determined 
from  this  known  bending  moment  by  the’ use  of  equation  (42).  If  the  mid¬ 
chords  of  the  two  airplanes  are  equal,  the  time  interval  may  be  taken 
equal  and  equation  (42)  reduces  to 

|K|2.|E][D]-"]2[[D]H-"|Ki]^  (43) 


SUMMARY  OF  CALCULATION  PROCEDURE 


As  a  convenience,  a  summary  of  the  basic  steps  necessary  for  calcu¬ 
lating  the  response  of  an  airplane  to  a  gust  is  given  as  follows: 

For  accelerations  and  displacements: 

(1)  With  the  use  of  the  fundamental  mode,  wing  plan  form,  mass 

distribution,  calculate  the  quantities  |Iq,  p.^,  X,,  r^,  and  as 

given  by  equations  (I8). 

(  . 

(2)  Choose  the  time  Interval  e.  A  convenient  rule  of  thumb  is 
€  ~  but  for  most  cases  6=1  should  give  satisfactory  results. 

(3)  Determine  values  of  the  unsteady-lift  function  0  =  1  _  $  at 
successive  multiple  intervals  of  e.  (See  fig.  1.)  Also  determine 
corresponding  values  of  the  gust- force  integral  f(s),  equation  (27b). 

As  an  aid,  curves  for  f(s)  are  presented  in  figure  1  for  the  sharp-edge 
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p 

gust  and  in  figure  2  for  various- length  sine  gusts,  sine  gusts,  and 
triangular  gusts.  (The  curves  in  fig.  1  have  been  obtained  from 
eqs.  (9)  and  (lO).  These  approximations,  although  rather  accurate  for 
the  lower  values  of  s,  are  noted  to  cross;  actually,  they  should  not 

cross  and  are  known  to  have  the  same  asymptotic  approach  to  unity. ) 

(4)  From  the  following  definitions: 


=  Mo 


set  up  ttie  following  matrices: 


(m  >  1) 


(m  >  1) 


(m  >1) 


Then  calculate  the  matrix 


[D]  =  i[A][c]  *  [b] 

^0 

(5)  Solve  for  the  values  of  P  (which  equals  from  equa¬ 

tion  (33) ^  by  the  method  outlined  following  equation  (33).  (See 

eq.  (3^).)  The  values  of  z-^  and  a  (which  equals  Zq")  can  then  be 

calculated  from  equations  (31)  and  (32). 

For  bending  moment; 

(6)  In  order  to  compute  bending  moment,  determine  r^,  and 

as  given  by  equations  (24),  where  M^,^,  and  in  these 

equations  depend  on  the  particular  wing  station  being  considered  and  are 
given  by  equations  (21). 
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(7)  Determine  bending  moment  by  use  of  equation  (23)  with  the 
values  of  response  already  established.  This  equation  may  be  applied 
directly  to  any  desired  time  value.  Maximum  bending  moment  usually 
occurs  very  close  to  the  time  when  is  a  maximum. 

EXAMPLES 


Example  A.-  In  order  to  provide  an  illustration  and  give  an  idea 
of  the  accuracy  of  the  present  analysis^  the  response  to  a  sharp-edge 
gust  of  the  two-engine-airplane  example  considered  in  reference  1  was 
determined.  The  weight  distribution  over  the  semispan^  the  wing-chord 
distribution^  and  the  fundamental  bending  mode  are  shown  in  figures  3^ 
ky  and  The  frequency  and  deflection  of  the  fundamental  mode  were 
calculated  by  the  method  given  in  reference  10.  The  solution  is  made 
for  a  forward  velocity  of  210  miles  per  hour  and  a  gust  velocity  of 
10  feet  per  second. 

The  lift-curve  slope  used  in  reference  1  was  5*^1;  to  be  consistent, 
the  same  value  was  used  here.  Furthermore,  the  unsteady-lift  function 
used  for  a  change  in  angle  of  attack  in  the  exan^le  presented  in  refer¬ 
ence  1  was  given  by  the  equation 


(1  -  =  1  -  0.36le"°*38ls 

rather  than  by  equation  (9).  Thus,  this  equation  was  also  used  here. 

The  gust  unsteady-lift  function  used  was  that  given  by  equation  (lO). 

The  various  physical  constants  and  the  basic  response  and  bending- 
moment  parameters  are  given  in  table  1;  the  values  of  the  unsteady-lift 
function  and  the  values  of  the  gust  force  are  listed  in  table  2.  The 

matrices  [A],  [B],  and  [c]  used  in  the  solution  are  given  in  table  3* 

The  solution  for  response  is  shown  in  figure  6(a)  where  the  deflec¬ 
tion  coefficients  Sq  and  a^  in  inches  are  plotted  against  distance 

traveled  in  half-chords.  The  corresponding  deflection  quantities  for 
the  example  given  in  reference  1  were  deteimined  and,  for  comparison, 
are  also  shown  in  the  figure.  A  similar  comparison  is  made  in  fig¬ 
ure  6(b)  for  bending  stresses  at  the  fuselage  and  engine  stations, 
stations  0  and  1  from  reference  1.  The  agreement  is  seen  to  be  good. 

Example  B.-  A  second  example  is  included  in  order  to  illustrate  one 
means  by  which  the  method  may  be  used  to  evaluate  the  influence  of  bending 
flexibility  upon  the  response  to  a  gust.  The  physical  characteristics  for 
the  airplane  considered  in  this  example  are  listed  also  in  table  1,  and 
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equation  (9)  is  used  for  the  function  1  -  $  instead  of  the  values 
given  in  table  2.  Maximum  values  of  the  bending  moment  that  develops 
at  the  fuselage  station  during  flights  through  sine  gusts  of  various 
lengths  have  been  determined,  both  for  the  aiiplane  considered  flexible 
and  for  the  aiipilane  considered  rigid.  The  results  are  shown  in  fig¬ 
ure  7(a)  where  maximum  values  of  the  bending-moment  parameter  K  are 
plotted  against  gust-gradient  distance  H.  The  difference  between  the 
two  curves  represents  the  increase  in  bending  moment  due  to  effects  of 
wing  bending  flexibility.  By  taking  the  ratio  of  K  for  the  flexible 
case  to  K  for  the  rigid  case,  a  type  of  dynamic  response  factor  is 
formulated  which  gives  a  direct  measure  of  the  influence  of  wing  flexi¬ 
bility.  This  ratio  is  designated  and  is  shown  in  figure  7(b).  As 

an  example  of  the  significance  of  this  plot,  the  value  of  =  I.16 

at  H  =  5  means  that  flexibility  results  in  a  l6-percent  dynamic  over¬ 
shoot  in  the  stress  from  the  value  that  would  be  obtained  at  H  =  5  on 
the  basis  of  a  rigid-body  analysis.  It  may  be  seen  also  that  the  value 
of  7j^  is  approximately  unity  for  values  of  H  =  10  and  greater;  there¬ 
fore,  in  this  range  of  gust-gradient  distances  a  rigid-body  treatment 
would  be  sufficient  for  this  aiip)lane. 


DISCUSSION 


The  derivation  presented  herein  is  intended  to  provide  a  convenient 
engineering  method  for  calculating  the  response  of  an  airplane  to  a  gust 
where  wing  bending  flexibility  is  included.  The  method  is  believed  to  be 
well-suited  for  making  trend  studies  which  evaluate,  for  example,  the 
effect  on  response  of  such  factors  as  mass  distribution,  speed,  and 
altitude.  Althou^  the  unsteady-lift  functions  for  two-dimensional 
unsteady  flow  are  presented,  the  method  is  general  enou^  so  that  the 
unsteady-lift  functions  for  finite  aspect  ratio,  for  subsonic  compress¬ 
ible  flow,  and  for  supersonic  flow  may  be  used  as  well.  (See  refs.  7 
and  11  to  15. ) 

Since  the  numerical  method  is  based  on  an  Integration  procedure,  it 
possesses  the  desirable  feature  that  a  fairly  large  time  interval  may  be 
used  and  good  accuracy  still  obtained.  As  an  accuracy  test,  solutions  of 
equations  (16)  and  (17)  were  made  for  several  cases  by  the  exact  Laplace 
transform  method  as  well  as  by  the  numerical  process,  in  which  process 
the  time  interval  was  selected  according  to  the  rule  of  thumb  suggested. 
When  the  results  were  plotted  to  three  figures,  the  difference  between 
the  two  solutions  was  barely  discernible. 

Additional  bending  modes  could  be  included  in  the  analysis  but  this 
refinement  is  really  not  warranted.  Some  calculations  made  with  addi¬ 
tional  modes  gave  results  which  differed  only  slightly  from  the  results 
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obtained  when  only  the  fundamental  mode  was  used.  The  good  agreement 
of  results  obtained  in  example  A  with  the  results  obtained  by  the  more 
precise  method  given  in  reference  1  also  illustrates  this  point. 
Furthermore^  if  additional  degrees  of  freedom  are  to  be  used,  it  would 
appear  more  important  to  include  wing  torsion  and  airplane  pitch.  The 
extent  to  which  torsion  influences  the  results  is  probably  governed  most 
by  the  nearness  to  the  flutter  speed.  The  importance  of  airplane  pitch 
is  probably  governed  most  by  the  gust  length;  some  investigations  dealing 
with  pitch  have  indicated  that  except  for  very  light  wing  loadings  the 
pitch  of  the  airplane  does  not  influence  the  results  appreciably  until 
gust  lengths  of  from  20  to  30  chords  or  larger  are  involved.  Thus,  the 
present  analysis,  although  limited  to  the  degrees  of  freedom  of  vertical 
motion  and  wing  bending,  should  probably  be  sufficiently  satisfactory 
for  speeds  near  the  cruising  speed  and  for  gust-gradient  distances  up  to 
approximately  10  chords. 

The  analysis  may  be  useful  in  assessing  the  significance  of  wing 
flexibility  in  the  technique  of  measuring  gust  intensity  by  means  of  an 
airplane.  In  this  technique  gust  severity  is  usually  measured  by  means 
of  an  accelerometer  placed  at  the  center  line  of  the  airplane.  In  order 
to  obtain  a  rough  idea  of  whether  flexibility  may  have  some  effect  on 
this  measurement,  calculations  for  the  maximum  accelerations  at  the 
center  line  and  for  the  maximum  acceleration  at  the  nodal  points  (the 
true  center- of- gravity  acceleration)  may  be  made  for  various  assumed 
gust  lengths.  A  comparison  of  these  computed  maximum  acceleration  values 
should  give  some  idea  as  to  the  extent  to  which  wing  flexibility  may 
alter  the  measurements  in  actual  flight. 


CONCLUDING  KEMAEKS 


The  analysis  presented  herein  for  the  response  of  an  airplane  to  a 
gust  should  provide  a  useful  means  for  evaluating  the  effects  of  wing 
flexibility.  A  convenient,  but  accurate,  numerical  solution  of  the 
response  equations  is  developed  which  is  well-suited  for  trend  studies 
such  as  the  evaluation  of  the  effects  of  mass  distribution,  speed, 
altitude,  and  similar  factors. 

As  indicated  by  an  example,  the  method  gives  good  agreement  with 
the  results  of  the  more  precise  but  more  lengthy  recurrence  matrix 
method  of  NACA  Rep.  1010. 

The  method  permits  the  evaluation  of  a  gust  causing  a  known  response. 
A  procedure  is  given  wherein  the  known  gust  response  of  one  airplane  may 


be  used  directly  to  determine  what  the  response  would  be  for  another 
airplane  flying  throu^  the  same  gust  condition. 


Langley  Aeronautical  Laboratory^ 

National  Advisory  Committee  for  Aeronautics^ 
Langley  Fields  Va. ,  May  21,  1952. 
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APPENDIX 


DERIVATION  OF  EQUATION  RELATING  DISPLACEMENT  TO 
PREVIOUS  SUCCESSIVE  VALUES  OF  ACCELERATION 

In  this  appendix^  a  derivation  is  given  of  equation  (31)  which  gives 
the  value  of  displacement  in  terms  of  successive  past-history  values  of 
acceleration.  Suppose  that  the  second  derivative  (acceleration)  of  a 
function  is  approximated  by  a  succession  of  straight-line  segments  as 
shown  in  the  following  sketch: 


where  the  segments  cover  equal  intervals  e  of  the  abscissa  s  and  the 
initial  condition  that  Zq"  =  0  is  assumed  to  apply.  If  a  dummy  origin 

is  now  considered  at  the  station  m  -  the  segment  between  sta¬ 
tions  m  -  1  and  m  may  be  represented  by  the  equation 


Two  successive  integrations  give  the  relations  for  and  Zjjj  as 

follows: 
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where  the  constants  of  integration  conditions 

for  the  interval)  have  been  introduced.  If  s  is  set  equal  to  e  in 
these  two  equations^  the  following  equations  result: 

"'m  =  +  ^'m-1  (“) 

2;^=^^2ijji  +  '^Z  +  ^V-1^  ^m-l 


From  these  two  equations  the  values  of  and  z^  at  any  time 

interval  may  be  given  in  terms  of  the  second  derivative  at  all  previous 
time  intervals.  For  example^  with  initial  conditions  of  z”q  =  z’q  =  0^ 

equation  (Al)  becomes  for  m  =  1 


(A3) 


and  for  m  =  2 


2*2  "  2 (^"2  ^"1)  ^'1 

Combining  this  equation  and  equation  (A3)  results  in  the  relation 

z'2  =  €(z”l  +  1  z"2) 


This  process  may  be  carried  through  for  each  of  the  time  stations  to 
yield  the  following  general  equation  for 


z 


m 


+  z 


+  z 


+  z 


m-1 


(A4) 


which^  of  course^  is  the  trapezoidal  approximation  of  the  area  under  the 
z" -curve.  Equation  (A2)  for  Zj^  may  be  treated  similarly  and  it  is 
found  that  the  general  equation  for  z^  may  be  written 


.2 


■^m 


(m  -  l)z^^2_  +  (m  -  2)z*’2 


2z” 
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■2+2  m-I  + 


1  2" 
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(A5) 


34 


NACA  TN  2763 


This  equation  thus  gives  the  displacement  at  any  time  station  in  terras 
of  the  accelerations  at  all  previous  time  stations. 

It  may  be  noted  that,  if  higher-order  segments  (parabolic  or  cubic) 
had  been  used  instead  of  straight-line  segments  to  approximate  the 
second  derivative,  equations  similar  in  form  to  equations  (a4)  and  (A5) 
would  also  result.  For  most  practical  purposes,  however,  the  accuracy 
of  equation  (A5)  is  sufficiently  good  as  long  as  the  interval  e  is 
chosen  so  that  the  straight-line  segments  roughly  approximate  the  second 
derivative. 
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TABLE  1.-  PHYSICAL  CHARACTERISTICS  FOR  EXAMPLE  AIRPLANES 


Example  A 

Example  B 

W,  It . 

37,450 

33,450 

S,  ft^  . 

870 

870 

b^  in . 

1120 

1120 

in . . 

154 

164 

p,  Ib/ft3 . . . 

0.0765 

0. 0765 

V,  ft'/sec . 

308 

374 

ft/sec  . 

10 

€  /sec  . 

0.0208 

0.0183 

^alf- chords . 

1.0 

1.0 

a . 

5-41 

6.28 

^^0 . 

64.16 

46.8 

. 

0.9045 

0.748 

\ . 

0.4353 

0.392 

^I . 

0.2181 

0.225 

"^2 . . . 

0.1358 

0.143 

P  J  fuselage  station . 

0.452 

0.457 

3  1  enerine  station . . 

. . 

0.  547 

1  fuselage  station . 

23.49 

15.94 

1 eneine  station  .  s  . 

10.19 

^  J  fuselage  station . 

3.665 

2.555 

[engine  station . 

3.  391 

*z  -3jf\iselage  station . 

0.00537 

0.00537 

^ngine  station . 

0.00669 

z  here  denotes  distance  from  neutral  axis  to  extreme  fiber. 
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TABLE  2.-  1  -  $  ORDINATES  AND  GUST-FORCE  ORDINATES 

FOR  SHARP-EDGE  GUST,  e  =  1.0 


m 

ejn  or  (1  - 

f  or 

0 

0.6390 

0 

1 

.7534 

.377 

2 

.8315 

.547 

3 

.881+9 

.635 

h 

.9211+ 

.692 

5 

.9463 

.736 

6 

.9633 

.771 

7 

.9749 

.798 

8 

.  9829 

.821 

9 

.9883 

.81+5 
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TABLE  3 . -  MATRICES  USED  IN  EXAMPLE  A 


A  Matrix 


64.799 

1.5068  64.799 

1.6630  1.5068 

1.7698  1.6630 

1.8428  1.7698 

1.8926  1.8428 
1.9266  1.8926 
1.9498  1.9266 
1.9658  1.9498 
1.9766  1.9658 


64. 799 

1.5068  64.799 

1.6630  1.5068 

1.7698  1.6630 

1.8428  1.7698 

1.8926  1.8428 
1.9266  1.8926 
1.9498  1.9266 


64. 799 

1.5068  64.799 

1.6630  1.5068 

1.7698  1.6630 

1.8428  1.7698 

1.8926  1.8428 


64. 799 

1.5068  64.799 
1.6630  1.5068 
1.7698  1.6630 


64. 799 

1.5068  64.799 


B  Matrix 


0.1394 
.  3286 
.3627 
.3860 
.4019 

0.1394 

.3286 

.3627 

.3860 

0.1394 

.3286 

.3627 

0.1394 

.3286 

0.1394 

.4128 

.4019 

.3860 

.3627 

.3286 

.4202 

.4128 

.4019 

.3860 

.3627 

.4252 

.4202 

.4128 

.4019 

.3860 

.4287 

.4252 

.4202 

.4128 

.4019 

.4311 

.4287 

.4252 

.4202 

.4128 

0.1394 

. 3286  0. 1394 

.3627  .3286  0.1394 

.3860  .3627  .3286  0.1394 

.4019  .3860  .3627  .3286  0.1394 


C  Matrix 


4.5367 

1.3954  4.5367 
2.2445  1.3954  4.5367 
3.0735  2.2445  1.3954 
3.8889  3.0735  2.2445 
^.6949  3.8889  3.0735 
5.4947  4.6949  3.8889 
6.2900  5.4947  4.6949 
7.0824  6.2900  5.4947 
7.8726  7.0824  6.2900 


4.5367 

1.3954  4.5367 
2.2445  1.3954  4.5367 
3.0735  2.2445  1.3954 
3.8889  3.0735  2.2445 
4.6949  3.8889  3.0735 
5.4947  4.6949  3.8889 


4.5367 

1.3954  4.5367 
2.2445  1.3954  4.5367 
3.0735  2.2445  1.3954  4.5367 


the  gust  force  f(s)=\/^(s). 
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Figure  2  r  Continued. 
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Figure  2.-  Continued. 
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Figure  2r  Concluded. 
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Figures.-  Semispan  weight  distribution  for  the  two-engine  airplane  of 

example  A. 


Figure 4.- Wing  chord  distribution  for  airplane  of  example  A. 


=  20.9  radians  per  second. 
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Figure  6,- Response  of  example  airplane  A  to  a  lO-foot-per-second  sharp-edge  gust. 

V=2I0  miles  per  hour. 
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(b)  Dynamic  response  factor  for  moment. 

Figure  7  r  Bending  moment  and  dynamic  response  factor  for  airplane 
of  example  B  due  to  flight  through  sine  gusts. 
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